Collisions of solitons and vortex rings in cylindrical Bose-Einstein condensates 
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Interactions of solitary waves in a cylindrically confined Bose-Einstein condensate are investigated 
by simulating their head-on collisions. Slow vortex rings and fast solitons are found to collide elas- 
tically contrary to the situation in the three-dimensional homogeneous Bose gas. Strongly inelastic 
collisions are absent for low density condensates but occur at higher densities for intermediate ve- 
locities. The scattering behaviour is rationalised by use of dispersion diagrams. During inelastic 
collisions, spherical shell-like structures of low density are formed and they eventually decay into 
depletion droplets with solitary wave features. The relation to similar shells observed in a recent 
experiment [Ginsberg et al. Phys Rev. Lett. 94, 040403 (2005)] is discussed. 

PACS numbers: 03.75.Lm, 47.32.Cc, 47.37.+q 



The observation of dark solitons and vortex rings in 
a series of experiments in Bose-Einstein condensates 
(BECs) [j], 0, EJ H has shown that this is an excellent 
physical system for the study of nonlinear waves. Theo- 
retical work has shown that the latter have the form of 
solitons, vortex rings, and solitonic vortices [1,0, ^ n 
the recent experiment of Ginsberg et al. collisions be- 
tween solitary waves were observed. While vortex rings 
and solitons were robust in many collision events, in some 
cases shell structures of low particle density were ob- 
served, which had not been predicted before. 

In the one-dimensional (ID) nonlinear Schrodinger 
equation (NLS), which describes the ID Bose gas, soli- 
ton collisions are elastic, that is no energy is radiated and 
the outgoing solitons are the same as the colliding ones 
[To| . In the homogeneous three-dimensional (3D) Bose 
gas, solitary waves have the form of vortex rings or rar- 
efaction pulses Their interactions are important for 
the understanding of superfluid turbulence. In contrast 
to the ID case, collisions are generally inelastic. Large 
vortex rings annihilate when they collide head-on by in- 
creasing their radius and radiating phonons 0, ^| . Col- 
lisions at oblique angles or with an impact parameter re- 
sult in vortex-line reconnections or produce Kelvin-wave 
radiation 0, ^| . Solitary wave collisions in tra ppe d 2D 
systems were previously considered in Refs. fl6l Il7j . 

A model system that connects the ID and 3D solitary 
waves is given by the cylindrically confined BEC. Study- 
ing elementary collision processes in this model is thus of 
fundamental theoretical interest and furthermore it leads 
to a deeper understanding of the experimentally observed 
structures. 

In this Letter we present a detailed numerical and theo- 
retical study of head-on collisions between solitary waves 
in a cylindrically confined BEC. We test the robustness 
of solitary waves when they interact. The detailed be- 
haviour of solitary waves under collision for low density 
condensates is found to be similar to soliton dynamics 
in the NLS This is, however, only a limiting case 
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FIG. 1: The energy fraction AE, defined in Q. which is radi- 
ated during collision of solitons and vortex rings of type I and 
II (opposite direction) as a function of their initial velocity, for 
three values of the coupling constant 7 = n\a. The velocity 
is normalised to the speed of sound v s = 0.95, 1.29, 1.61 (in 
units of h/(ml p )) for 7 = 1,3,7, respectively. The symbols 
correspond to the simulations performed. 



for the present system, which also shows more compli- 
cated dynamical behaviour. For higher particle densities 
we find that the collision dynamics of solitary waves, in- 
cluding vortex rings, is very different than in the homo- 
geneous 3D Bose gas. As shown in Fig.^ we find elastic 
collisions for small and large velocities. At intermedi- 
ate energies, however, inelastic collisions often produce 
temporary spherical shells reminiscent of those observed 
experimentally |9|. In the following, we introduce the 
model and discuss the various regimes in detail. 

We assume a cylindrical and infinitely-elongated trap 
as in Ref. 7] with symmetry axis z. The Gross-Pitaevskii 
equation can be written in the dimensionless form 



~dt 
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FIG. 2: The excitation energy in units of nih 2 /(ml p ) ver- 
sus impulse Q in units of mh for axially symmetric solitary 
waves 7]. At 7 = 1 in panel (a) there is one family of solitons. 
Panel (b) for 7 = 7 shows two fundamental branches (solid 
lines) of vortex rings related by @. The vortex ring radius 
has a minimum at the cusps, where \v\ = 0.3v s . The im- 
pulse transfer during bounce-off collisions is indicated by the 
dashed arrow connecting initial solitary waves (filled symbols, 
see text) with collision products (open symbols). The vertical 
dotted lines indicate the energy loss. 



where p = ^/ x 2 + y 2 is the radial coordinate. The dimen- 
sionless coupling constant 7 = rt\a is the only parameter 
entering the equation, where a is the scattering length 
and n\ is the linear particle density along the symme- 
try axis z. Length is measured in units of the oscillator 
length l p = yjfi/moup, where m is the atomic mass and 
Up is the transverse trapping frequency. The unit of time 
is l/(jp. At z — > ±00 the wave function approaches the 
ground state in the transverse plane with jdz = and 
f\V\ 2 d 2 x = 1. 

The wave functions of solitary waves in the cylindrical 
trap are of the form 

9(x, y, z, t- v) = ^{x, y,z- vt) e~^\ (2) 

where v is the velocity of propagation and [i the chemical 
potential. For 7 < 1.5 only one family of solitary waves 
exists with velocities — v s < v < v s . They are related 
to the dark soliton of the NLS. We show the energy- 
momentum dispersion for 7 = 1 in Fig.[2K 7]. A modified 
momentum Q, called the impulse, has been defined so 
that the velocity is given by the slope of the curve v = 
dE/dQ Solitons with opposite velocities are obtained 
by the symmetry transformation 

* **, v -v, (3) 

where the star denotes complex conjugation. 

We numerically simulate the elementary process of a 
head-on collision of two solitons at 7 = 1 with wave func- 
tions related by Eq. We assume axial symmetry 
throughout the collision process and we typically use a 




FIG. 3: Plots of the particle density |^| 2 for solitary wave 
collisions on a scale from white at zero to black at maximum 
density for each plot. Each row shows five snapshots of a 
solitary wave collision with values of 7 = 1 (a,b), and 7 = 3 
(c), and velocities v as indicated. The time for each snapshot 
is shown in the figure. The spatial dimensions of the plots are 
8 x 12 for (a,b) and 8 x 16 for (c). The initial separation (at 
time t = 0) of the solitary waves is d = 12 in units of l p for 
all cases. 

lattice size 6 x 60 in the (p, z) plane and a lattice spac- 
ing 0.1 for both variables. We integrate Eq. in time 
using a Runge-Kutta method for various values of veloc- 
ity v of the initial solitons each calculated as in Ref. 
The outgoing solitons after collision are in all cases very 
similar to the incoming ones and no visible radiation is 
produced, thus collisions are elastic. A measure of the 
energy transformed to radiation during collision is 



where E{ is the energy of the incoming solitons and Ef 
is the energy of the outgoing ones. The energy Ef of 
the outgoing solitons is inferred by measuring their ve- 
locity. As seen in Fig.^ reaches a maximum 5% for 
intermediate velocities v « 0.5v s (v s = 0.95). 

In Fig. we present five snapshots of the simulation 
for solitons with initial velocity v = 0.2(0.21v s ). The 
solitons decelerate as they approach and interact, they 
reach a minimum separation and finally bounce back. 
This behaviour is typical for small velocities. One of the 
initial solitons is denoted by a filled circle in Fig.EK? and 
it appears that it effectively moves along the dispersion 
curve due to the collision process to the point denoted 
by an open circle at almost the same energy level. This 



picture gives a precise result for the final outcome of the 
collision, and it also describes gross features of the colli- 
sion process. It cannot give the detailed features of the 
process because the dispersion pertains only to isolated 
solitary waves. As the velocity of the initial solitons in- 
creases, they pass through each other during collision as 
is seen in Fig. [2b for v = 0.6(0.63v s ). In the dispersion 
diagram of Fig.EK, one of the solitons is represented by a 
filled square before and by an open square after collision. 
The energy difference between the two points is small, 
thus the collision is almost elastic. The overall picture 
at low 7 closely resembles the elastic soliton collisions in 
the integrable NLS 0. 

We repeat the simulations for a denser condensate by 
setting 7 = 3 in Eq. Vortex rings, albeit with a very 
inconspicuous ring structure, form now part of the soli- 
tary wave family. The energy radiated during collision 
reaches now a maximum 23% at v = 0Av s (v s — 1.29) as 
shown in Fig. Nevertheless, collisions for slow and for 
fast solitary waves are also in this case elastic and their 
behaviour resembles that for 7 = 1 0. In Fig. Eh we 
show five snapshots of the simulation for v = 0.4(0.31^). 
The second snapshot shows the formation of fully fledged 
vortex rings at the time of collision which do not exist 
as isolated solitary waves. In fact we have noticed the 
transient formation of vorticity during solitary- wave col- 
lisions in many cases in this work. Similar observations 
were also made in the simulations of the Ginsberg experi- 
ment 9] and for simulated collisions of rarefaction pulses 
with vortex lines in the homogeneous BEC The 
outgoing waves are seen to be of the soliton type in the 
fourth snapshot. However, these are later spontaneously 
transformed to depletion droplets, i.e. strongly localised 
regions of low density, as seen in the last snapshot. At 
later times not shown in the figure, the droplets trans- 
form into soliton-like structures and periodically revive 
two to three times until we can no longer follow the sim- 
ulation due to lattice limitations. The speed of the wave 
is almost constant at 0A8v s . 

We further study how the picture changes as the cou- 
pling increases. The energy-momentum dispersion for 
axially symmetric solitary waves changes substantially 
for 7 > 4, and it is shown for 7 = 7 in Fig. ^jp |g. It 
has two fundamental branches (solid lines in the figure) 
which are related by Eq. and they both contain axi- 
ally symmetric vortex rings with velocities of both signs. 
For example, for the branch originating on the left we 
have 1.61 = v s > v > —v\ = — 0.30v s where —v\ is the 
velocity at the cusp. For velocities close to v s , the soli- 
tary waves have the structure of gray solitons. For lower 
velocities they are vortex rings whose radius decreases as 
v decreases (as v — > —v\). In the following we need to 
distinguish between the vortex rings on the left and those 
on the right of the dispersion maximum. We call the vor- 
tex rings with v > type I, and those with —v\ < v < 
type II. Finally, the symmetry relation © gives a second 




FIG. 4: Particle density |^| 2 for solitary wave collisions for 
7 = 7. The frames are 8 x 12. 



branch in Fig. [2)3 with mirror imaged type I and II rings. 
The branch denoted by a dotted line in Fig. (2t> contains 
unstable solitary waves of higher energy and will not be 
discussed further here. 

We simulate the collision of two counter-propagating 
vortex rings of opposite circulation related by Eq. (0. 
We first discretized the wave function in a 3D lattice in 
Cartesian coordinates. After checking in test runs that 
the axial symmetry was not broken during collision we 
proceeded to extensive simulations assuming axial sym- 
metry throughout the process. The energy radiated due 
to collision for a coupling constant 7 = 7 is shown in 
Fig. H The maximum is 60% and occurs for initial ve- 
locity v = 0.37v s , (v s = 1.61). On the other hand, colli- 
sions at low and high velocities appear to be elastic, i.e., 
the outgoing solitary waves have almost the same energy 
as the incoming ones. The three curves in Fig. for 
7 = 1,3, and 7 show that collisions are elastic for all cou- 
plings when the colliding solitary waves are slow or when 
their velocity is close to the speed of sound [l£j- For in- 
termediate velocities energy is radiated during collision 
and this is higher as the coupling increases. 

Fig-Bt shows the collision of vortex rings with an initial 
velocity v = 0.3(0. 19v s ). The initial rings are denoted by 
a filled circle in the dispersion (Fig. 03 ) and they are type 
I. They reach a minimum distance, their radius decreases, 
and subsequently they bounce back. The outgoing rings 
of type II still have the original circulation and are de- 
noted by an open circle on the dispersion curve. The 
change in nature of the rings is apparent in the change 
of the ring radius between the first and last snapshots in 
Fig. The energy radiated during this process is 3%. It 
appears that each ring has effectively moved along one of 
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the branches of the dispersion curve. This picture gives 
a precise result for the final outcome of the collision, but 
it also gives a faithful representation of the gross features 
throughout the collision process. The almost elastic col- 
lisions of slow vortex rings seen here are in stark contrast 
to collisions in the 3D bulk where vortex rings increase 
their radius and eventually annihilate 12]. They are in- 
hibited to do this in our model due to the transverse 
confinement. 

We may choose the initial vortex rings to be type II 
but still have initial velocity \v\ = 0.3(0.19^) (Fig.gJ)). 
In this case the collision causes them to increase their ra- 
dius, and they are eventually transformed to type I rings 
moving in opposite directions. The process is essentially 
the reverse of the one in Fig. 

The arguments pertaining to the dispersion curve can 
be employed to rationalise the observed behaviour for 
higher velocities (lower energies). One can thus imagine a 
collision where the initial type I vortex rings have energy 
slightly lower than the energy at the cusp, as for example 
at the point denoted by a filled square in Fig. [3d (corre- 
sponding to v = 0.5(0.31v s )). The bounce-back process 
of vortex rings cannot occur in this case. Instead, inter- 
mediate structures are formed, as discussed below, and 
the collisions become highly inelastic. In the example of 
Fig. 03 the final waves, denoted by an open square, have 
substantially lower energy than the initial rings. This 
mechanism of inelastic collision suggests that the maxi- 
mum energy loss occurs at intermediate velocities. 

A further example of a collision for an initial veloc- 
ity 0.55 (0.34v s ) for type I rings is shown in Fig. QJ. 
The vortex rings form a shell-like object of low density 
as they collide, as seen in the third entry of the figure. 
Since we have assumed axial symmetry, the shell is actu- 
ally almost spherical, and strongly reminiscent of those 
observed in Ref. 9]. However, the shells reported here 
are produced by an elementary head-on collision process 
which is significantly simpler than the process reported 
in the experiment where many nonlinear wave structures 
interact simultaneously while the condensate is expand- 
ing. The resemblance between the observed structures 
suggests that inelastic collisions of axisymmetric nonlin- 
ear waves is the fundamental process underlying the gen- 
eration of spherical waves observed experimentally 9] 

As the rings move away to opposite directions they 
form depletion droplets (similar to those seen in Fig. 
which are shown in the last entry of the figure. These 
travel coherently for a distance of approximately 10 units 
and they eventually seem to decay into solitons. The 
depletion droplets are distinctly different than solitary 
waves theoretically studied in a confined BEC but they 
are reminiscent of the rarefaction pulses in the homoge- 
neous Bose gas 0, [H| • 

The collision in Fig. is highly inelastic in the sense 
that the main outgoing waves carry only part of the to- 
tal energy while the rest of the energy is radiated away. 



However, it is possible that some of the remaining energy 
is carried by shallow gray solitons. 

At small energies (v close to v s ) the solitary waves have 
a gray-soliton character. They pass through each other 
during collision much like the situation for small cou- 
plings. In Fig.Eb the solitary wave with v = 1.0(0.62i> s ) 
is denoted by a filled triangle. After the collision it has 
only slightly lower energy (open triangle). Thus, colli- 
sions for large velocities are elastic as is also shown in 
Fig. |H 1$. 

Concluding, we have shown that elastic collisions of 
solitary waves can occur in cylindrical BECs. Shell struc- 
tures reminiscent of recent experimental observations 
were shown to arise already in elementary inelastic head- 
on collisions. Possible extensions of the present work 
to larger densities, non- axisymmetric solitary waves, and 
beyond the head-on case will provide further valuable in- 
sights into the dynamics of nonlinear waves in confined 
BECs. 
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